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ABSTRACT 

Biogenesis of trans-acting siRNAs (tasiRNAs) is 
initiated by miRNA-directed cleavage of TAS gene 
transcripts and requires RNA-dependent RNA poly- 
merase 6 (RDR6) and Dicer-like 4 (DCL4). Here, we 
show that following miR173 cleavage the entire 
polyadenylated parts of Ambidopsis TAS1a/b/c and 
TAS2 transcripts are converted by RDR6 to 
double-stranded (ds)RNAs. Additionally, shorter 
dsRNAs are produced following a second cleavage 
directed by a TASfc-derived siRNA. This tasiRNA 
and miR173 guide Argonaute 1 complexes to 
excise the segments from TAS2 and three TAS1 
transcripts including TAS1c itself to be converted 
to dsRNAs, which restricts siRNA production to a 
region between the two cleavage sites. TAS1c is 
also feedback regulated by a c/'s-acting siRNA. We 
conclude that TAS1c generates a master siRNA that 
controls a complex network of TAS1/TAS2 siRNA 
biogenesis and gene regulation. TAS1/TAS2 short 
dsRNAs produced in this network are processed 
by DCL4 from both ends in distinct registers, 
which increases repertoires of tasiRNAs. 



INTRODUCTION 

Micro RNAs (miRNAs) regulate gene expression by 
silencing their target genes through mRNA cleavage 
and/or translational repression (1,2). In plants, some 
miRNA targets spawn 21-nt secondary siRNAs (3,4), a 
subset of which, called trans-acting siRNAs (tasiRNAs), 
can also target genes in trans. The tasiRNA-targeted genes 
can in turn spawn tasiRNAs (5), thus leading to silencing 
cascades which are predicted to form complex networks 



with distinct regulatory functions (6). miRNAs, tasiRNAs 
and secondary siRNAs produced from natural antisense 
transcripts are involved in regulation of plant develop- 
ment, physiology and stress responses (7-9). 

The biogenesis of tasiRNAs requires RNA-dependent 
RNA polymerase 6 (RDR6) which converts a portion of 
target RNA upstream or downstream of the miRNA/ 
siRNA cleavage site to double-stranded RNA (dsRNA). 
The resulting dsRNA is then processed by Dicer-like 4 
(DCL4) from the miRNA-cleavage end into 21-nt 
in-phase siRNA duplexes with 2-nt 3'-overhangs (4,10,1 1). 

The Arabidopsis thaliana genome contains four families 
of tasiRNA-generating TAS genes, classified based on 
miRNA-dependence, sequence similarity and target 
gene identity. miR173-dependent TASla/b/c and TAS2 
tasiRNAs and miR828-dependent TAS4 tasiRNAs are 
encoded downstream of the miRNA cleavage site and 
target the genes of pentatricopeptide repeat and MYB 
transcription factor families, respectively (3,4). miR390- 
dependent TAS3a/b/c tasiRNAs are encoded upstream 
of the miRNA-cleavage site and target the genes of 
auxin response factor family (4,10). 

One of the open questions concerning secondary siRNA 
biogenesis is which factors route some (but not all) 
miRNA/tasiRNA target transcripts into the RDR6/ 
DCL4 pathway. For the Arabidopsis TAS3 family, a 
second, non-cleavable, miR390-binding site located 
upstream of the miR390 cleavage site is essential for 
tasiRNA biogenesis (12) and the region of TAS3a/b/c 
transcripts between the two miR390 sites spawns most 
of secondary siRNAs including tasiRNA species (4). 
Interestingly, TAS3 transcripts from moss and pine are 
cleaved at both miR390-target sites (12). These results, 
together with a finding that secondary siRNAs more 
often emanate from the genes with two known miRNA/ 
tasiRNA target sites, suggested a two-hit trigger model for 
secondary siRNA biogenesis (12,13). However, this model 
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left unexplained a robust production of siRNAs from 
miR173-dependent TAS genes. Efforts to identify a 
second hit by a known small RNA or some other tran- 
script feature required for tasiRNA production have failed 
so far (14). 

miRNAs and siRNAs exert their silencing functions in 
association with Argonaute (AGO) proteins that sort 
sRNAs based on size and 5'-terminal nucleotide 
(13,15,16). Thus, AGOl binds mostly 21-nt RNAs with 
5'-terminal uridine (5'U), while AG02 binds mostly 21-nt 
RNAs with 5'-terminal adenine (5' A). Interestingly, both 
miR173 and miR390 that initiate tasiRNA biogenesis 
deviate from this general rule. miR390 is a 21-nt 5'A- 
RNA but it is specifically associated with AG07 that 
plays an essential role at both cleavable and non-cleavable 
TAS3 target sites (13,17). miR173 is a 22-nt 5'U-RNA 
that associates with AGOl to guide the cleavage of 
TASla/b/c and TAS2 transcripts (14). The 22-nt size of 
miRNA is sufficient to trigger RDR6-dependent siRNA 
production from miRNA target transcripts including 
miR173 targets (18,19). A more recent study, however, 
indicates that asymmetrically positioned bulged bases in 
a miRNA:miRNA* duplex are more important than the 
miRNA length for the initiation of RDR6-dependent 
siRNA biogenesis (20). 

Here we used blot hybridization and circularization- 
reverse transcription-PCR (cRT-PCR) to characterize 
the structures of dsRNAs originating from Arabidopsis 
miR173-dependent TAS genes and immunoprecipitation 
(IP) to identify tasiRNAs associated with AGO 
complexes. Our study uncovers a complex sRNA 
network connecting TAS2 and TAS1 genes which is 
controlled by miR173 and a TASlc-derived tasiRNA. 



MATERIALS AND METHODS 

Arabidopsis wild-type, rdr6-15 and del mutant lines used in 
this study, their growth conditions and infection with 
CaMV, and the detailed protocols for total RNA prepar- 
ation and high-resolution blot hybridization have been 
described earlier (21,22). agol-25 and ago2-l were kindly 
provided by Dr Baulcombe. Seeds of taslc-1 
(SAIL_242_C04) were obtained from the Arabidopsis 
Biological Resource Center and the homozygous individ- 
uals were identified by genotyping. Probes used for 
RNA blot hybridization and genotyping are listed in 
Experimental Procedures in Supplementary Data. 
cRT-PCR mapping of TAS RNAs was performed as 
described earlier [(23); for details, see Methods in 
Supplementary Data]. 

The IP experiment with native AGOl -antibodies 
(Agrisera) and blot hybridization analysis of 
AGOl -associated sRNAs were performed as we described 
recently (24). The sRNA probes are listed in Methods in 
Supplementary Data. 

Transient expression experiments using Nicotiana 
benthamiana were performed as described previously 
(10). The 35S:TASlc construct was mutagenized by 
replacing the region between Pfo I and Pac I sites with 
synthetic gene fragments carrying the corresponding 



mutations. For further details, see Methods in 
Supplementary Data. 

RESULTS 

Long and short co-terminal dsRNAs originate from all 
miR173-targeted TAS genes 

Previously we found that infection with Cauliflower mosaic 
virus (CaMV) leads to stabilization of RDR6-dependent 
precursors of tasiRNAs in Arabidopsis (21). Using 
high-resolution gel separation of total RNA followed by 
blot hybridization with DNA oligonucleotide probes 
specific for TAS1 and TAS2 loci, we further characterized 
the stabilized sense (S) and antisense (AS) precursors 
(Figure 1). To investigate DCL requirements for dsRNA 
biogenesis, we compared CaMV-infected wild-type 
Arabidopsis (Col-0 and La-er) with double (dcl2/3 and 
dcl3/4) and triple (dcl2/3/4) Z)CL-mutant fines [all in 
Col-0 background; (22)]. All the mutants still express 
DCL1 that produces miRNAs. 

Given that tasiRNA siR255 is encoded in sense strands 
of all three TAS1 genes (la, lb and 7c) the probe comple- 
mentary to siR255 would detect their primary Pol II tran- 
scripts (pri-RNAs) and the sense strands of 
RDR6-dependent dsRNAs. The probe complementary 
to the star strand of siR255/siR255* duplex would 
detect only the antisense strand of RDR6-dependent 
dsRNAs. The two probes gave identical hybridization 
patterns in the range of 1 00-600 nt (Figure 1A), indicating 
that the RNAs detected within this range can potentially 
be sense and antisense strands of distinct dsRNA precur- 
sors of siR255/siR255*. To further classify these RNAs, 
we used the probes specific to the 21-nt sequences adjacent 
to the miR173 cleavage site (identical in both TASla and 
TASlb but not TASlc) and the probe specific to TASlb 
(Figure 1A). The resulting hybridization patterns allowed 
to identify the pri-RNAs and putative dsRNAs produced 
from each of the three TAS1 loci (Figure IB). Likewise, 
probes specific for the TAS2 sequences identified one 
major and several minor sense/antisense RNA pairs 
(Figure 1). 

Comparison of the hybridization patterns in wild-type 
versus <ic/-mutant lines indicated that most of the TAS1 
and TAS2 putative dsRNAs, except for the longest 
species, require DCL4 activity for their biogenesis. The 
longest sense/antisense pair derived from each locus 
over-accumulated in dcl2/3/4 (Figure 1A). The longest 
and shorter pairs from each locus, except for one TASlb 
and two TAS2 non-abundant pairs (Figure IB; designated 
DCL4 processing intermediates), all contain sequences 
adjacent to the miR173 cleavage site. Thus, most of the 
shorter dsRNAs are not intermediates of DCL4 process- 
ing of the longest dsRNAs from the miR173 end. The 
results described below are consistent with this conclusion. 

Confirming our previous results (22), CaMV-infected 
rdr6-15 mutant plants, which lack RDR6 activity, did 
not accumulate any TAS1 / TAS2-dehved sense/antisense 
RNA pairs. Only sense strands of the longest dsRNAs, 
i.e. miR173-cleavage products of TAS1 /TAS2 pri-RNAs 
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Figure 1. Blot hybridization analysis of dsRNA precursors of TASlajb/c and TAS2 tasiRNAs. (A) The panel shows blot hybridization analysis (5% 
PAGE) of total RNA from CaMV-infected (+) or mock-inoculated (-) wild-type (Col-0 and La-er), double (dcl2/3 and dd3j4) and triple {dcl2/3/4) 
DCL-deficient mutants. The scans represent the same blot successively hybridized with probes specific for TASIa/blc-, TASla/b-, TASlb-, TASlc- 
and 7^4S2-derived 21-nt siRNAs and U6 RNA as a loading control. Single-stranded RNA markers are shown on the left. The scheme shows location 
and polarity of the probed siRNAs (colored boxes with their names indicated) with respect to the miR173 cleavage site (indicated by arrow) on 
hypothetical dsRNAs from the four TAS genes. Identical siRNAs/probes, i.e. siR255/siR255* duplex encoded by TAS la, lb and lc and 3'D1(— ) 
siRNA of TAS la and lb, are indicated in pink and purple, respectively. Cyan, blue and orange indicate siRNAs/probes specific for TASlb, TASlc 
and TAS2, respectively. (B) Interpretation of the analysis shown in panel A. The lanes corresponding to Col-0 and dcl2/3/4 samples were cropped 
from the above blots of sense (S) polarity. Positions of pri-RNAs and dsRNAs are indicated by arrows for each TAS gene (la, lb, lc and 2). 
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(see below), could be detected by blot hybridization in 
rdr6-15 (Supplementary Figure SI). 

Polyadenylated 3' -products of miR173-cleaved TAS 
pri-RNAs are entirely converted by RDR6 to dsRNAs 

To precisely map TAS dsRNAs we used cRT-PCR that 
allows simultaneous sequencing of 5'- and 3'-termini of 
a given RNA [(23), Supplementary Figure S2A]. 
Pyrophosphatase treatment of total RNA, prior to RNA 
ligation, allowed circularization of transcripts not only 
with a 5'-monophosphate present on 3'-fragments of 
miRNA-cleaved TAS pri-RNA, but also with a 5'-triphos- 
phate likely present on an RDR6-synthesized complemen- 
tary RNA. cRT-PCR primers proximal to the termini 
of a given sense or antisense RNA (Figures 2 and 3; 
Supplementary Figure S2) were designed based on the 
sizes of putative dsRNAs estimated by blot hybridization 
(Figure IB). cRT-PCR products were gel-purified, cloned 
and sequenced (~ 10-30 individual clones for each). 

For each TAS, we first mapped the longest sense/anti- 
sense RNA pair from wild-type (Col-0) and dcl2/3/4 
plants. Typically, no abundant products could be 
obtained using total RNA from mock-inoculated 
Arabidopsis, whereas the CaMV-infected material-yielded 
abundant cRT-PCR products of expected size 
(Supplementary Figure S2). Along with our blot hybrid- 
ization (Figure 1), this indicates that tasiRNA precursors 
are efficiently processed by DCLs and do not accumulate 
at high levels, unless they are stabilized [in our case by the 
CaMV transactivator protein; (21)]. cRT-PCR analysis 
revealed that most of the TASla longest sense RNAs 
are the entire 3'-products of miR173-directed cleavage of 
the pri-RNA: they start with U379 (numbered from the 
cap site) and terminate with a poly(A) tail of up to 24 nt 
attached at one major site (A933) and several minor sites 
(Figure 2 and Supplementary Figure S2A). One to three 
terminal uridines were attached to some poly(A) tails. 
Most of the TASla longest antisense RNAs (RDR6 
products) terminate with A379 and a large fraction of 
those start with a poly(U) stretch of up to 26 nt attached 
to U933 or at nearby positions. Some antisense RNAs 
have short A-rich stretches at the very 5'-end 
(Supplementary Figure S2A). Thus, the antisense RNAs 
are essentially complements of the sense RNAs. These 
results demonstrate that RDR6 copies the entire 
3'-product of miR173-cleaved TASla pri-RNA into a 
complementary RNA. 

cRT-PCR analysis of the longest dsRNAs from TASlb, 
TASlc and TAS2 genes (Supplementary Figure S2B, S2C 
and S2D, respectively) led us to similar conclusions with 
some peculiarities in each case. For TASlb, two major 
(U753 and A822) and several minor poly(A) sites were 
identified (Figure 2 and Supplementary Figure S2B). The 
longest antisense RNAs from TASlb have very short 
U-stretches, UA-rich stretches or just a single uridine at 
the 5' termini, suggesting that long poly(A) tails are almost 
entirely removed from the 3'-cleavage product before its 
conversion to dsRNAs. For TASlc, two major poly(A) 
sites were mapped (A810 and G816), one of which 
(A810) was reported earlier (14). TASlc pri-RNAs of 



both types are used as templates for RDR6 following 
miR173-cleavage, since the counterpart antisense RNAs 
with long U-stretches at the 5'-end were found (Figure 2 
and Supplementary Figure S2C). The longest TAS2 
dsRNAs were composed of miR173-cleaved pri-RNA 
terminating at several minor and one major (A1031) 
poly(A) sites and their antisense counterparts (Figure 2 
and Supplementary Figure S2D). Notably, in this 
case the longest dsRNAs could be detected by cRT- 
PCR only in dcl2/3/4 but not in wild-type plants, 
which is consistent with the blot hybridization analysis 
(Figure 1). 



A second cleavage of miR173-targeted TAS pri-RNAs 
generates an additional template for RDR6 

cRT-PCR results summarized in Figure 3 show that sense 
RNA of the most abundant sense/antisense pair derived 
from TAS2 in wild-type plants starts with U418 and ends 
predominantly with G734. Thus, this RNA is a derivative 
of the miR173-cleaved pri-RNA 3'-product truncated by a 
second, downstream cleavage. The corresponding anti- 
sense RNA starts with A732 and ends with A418 
(Figure 3). Thus, RDR6 synthesis of a complementary 
RNA begins at the third nucleotide of the sense 
template, which creates a 2-nt 3'-overhang at the end of 
the resulting dsRNA, and stops at the template's last nu- 
cleotide. A large fraction of the antisense RNAs contained 
one to three non-template nucleotides, predominantly UU 
or UA (Supplementary Figure S2D). Since a tomato RDR 
activity was reported to add one non-template nucleotide 
to a fraction of the run-off antisense transcripts in vitro 
(25), we assume that the Arabidopsis RDR6 has a similar 
property in vivo. cRT-PCR analysis of the shorter TAS2 
dsRNAs in dcl2/3/4 plants which lack DCL4 activity 
revealed heterogeneous populations of sense and antisense 
RNAs with the sense RNAs starting at the miR173 
cleavage site but ending at various positions and the anti- 
sense RNAs mirroring their sense counterparts 
(Supplementary Figure S2D). We conclude that TAS2 
pri-RNA is cleaved twice, first at the miR173 site and 
second at the downstream site by a cleavage indirectly 
dependent on DCL4. The 317-nt RNA fragment 
between the two cleavage sites is converted by RDR6 to 
a 315-bp dsRNA (Figure 1). This dsRNA has a 2-nt 
3'-overhang at the second cleavage end and an up to 
3-nt 3'-overhang at the miR173 end. 

cRT-PCR analysis of TASlb short sense RNA also 
revealed a second, DCL4-dependent, cleavage site down- 
stream of the miR173 site: all the cloned RNAs start with 
A374 and end with G601 (Figure 3 and Supplementary 
Figure S2B). The corresponding antisense RNA starts 
either with C600 or U599 and ends predominantly with 
U374 followed in 50% cases by one or two non-template 
uridines (Supplementary Figure S2B). Thus, the most 
abundant precursor of TASlb tasiRNAs is 227- or 
226-bp dsRNA with a 2-nt 3'-overhang at the second 
cleavage end and an up to 2-nt 3'-overhang at the 
miR173 end. 
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Figure 2. The entire, polyadenylated 3'-products of miR173-cleaved pri-RNAs of TASlci, TASlb, TASlc and TAS2 genes are converted by RDR6 
to dsRNAs. Diagrammatic representation of TAS genes transcribed by Pol II into 5'-capped and 3'-polyadenylated pri-RNAs which are then cleaved 
by miR173-AG01 complexes. The 3'-products are converted by RDR6 to long dsRNAs. The regions of the miR173 cleavage site and the poly(A) 
sites of each TAS pri-RNA are enlarged. Arrows above the upper strand indicate positions of the 5'-terminus (miR173 site) and the 3'-terminus of at 
least one cloned cRT-PCR product representing the sense strands of TAS dsRNAs; the 3'-termini to which poly(A) is added are indicated in pink. 
Arrows below the complementary strand indicate the termini of one or more cloned cRT-PCR product representing the corresponding antisense 
strands of dsRNAs. Termini of the antisense RNA that begin with a poly(U) stretch are indicated in cyan. Thick arrows indicate the termini of the 
major sense and antisense RNAs; the number of specific cRT-PCR products/total products sequenced for each strand are indicated. Positions of 
cRT-PCR primers used for long dsRNA mapping are indicated above and below of the TAS gene body. 



The TASla sense/antisense pair of ~320bp is more 
abundant than the TASla longest dsRNA but less 
abundant than two shorter pairs derived from this locus 
(Figure IB). cRT-PCR mapping showed that this dsRNA 
is produced by the two-cleavage mechanism involving 
miR173-directed and DCL4-dependent cleavages. These 
cleavages excise a 318-nt fragment from TASla 
pri-RNA (from U379 to U696) which is copied by 
RDR6 to the complementary RNA starting either with 
U694 or C695 and ending with A379, often followed by 
one or two non-template nucleotides (Figure 3 and 
Supplementary Figure S2A). Thus, this dsRNA resembles 
the most-abundant TASlb dsRNA. 



TASlc-deri\ed tasiRNA directs a second cleavage of 
TAS2, TASla and TASlb pri-RNAs 

Our inspection of the second cleavage sites mapped on 
TAS2, TASlb and TASla pri-RNAs revealed a remark- 
able similarity: all these sites are near-perfectly comple- 
mentary to the r^S/c-derived siRNA 3'D6(-); this 
siRNA is encoded at the sixth processing cycle from the 
miR173-cleavage site in the antisense orientation (4,10). In 
all three cases, the cleavage occurs within a target 
pri-RNA/tasiRNA duplex between positions 10 and 11 
from the tasiRNA 5'-end (Figure 3), thus representing a 
conserved, AGO-mediated cleavage. Inspection of sRNA 
deep-sequencing data [http://asrp.cgrb.oregonstate.edu; 
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Figure 3. Two-cleavage mechanism for biogenesis of TASlci-, TASlb-, TASlc- and 7M,S2-derived dsRNAs. Diagrammatic representation of TAS 
genes transcribed by Pol II into pri-RNAs which are then cleaved at the two sites by miR173- and TASlc 3'D6(— ) siRNA-AGOl complexes, 
respectively. The excised fragment between the two-cleavage sites is converted by RDR6 to dsRNA. The cleavage sites for each pri-RNA and the 
siRNA/target pri-RNA duplexes are shown in the expanded regions. Additional complementary sites for TASlc 3'D6(— ) siRNA in the intervening 
regions of TAS2 and TASlci pri-RNAs as well as the origin site of TASlc 3'D6(— ) are also shown. Colored arrowheads indicate positions of the 
5'-terminus (miR173 site; in blue) and the 3'-terminus [TASlc 3'D6(— ) site; in red] of cloned cRT-PCR products representing the sense strands of 
TAS dsRNAs. The 5'- and 3'-terminal sequences of cRT-PCR clones representing the corresponding antisense strands of TAS dsRNAs are shown 
below the duplexes and highlighted in cyan. The number of specific cRT-PCR products/total products sequenced for each strand is indicated. 
Positions of potential cleavage sites for 3'D6(— ) siRNA are indicated by open arrowheads. Positions of the cRT-PCR primers used for the short 
dsRNA mapping are indicated above and below of the TAS gene body. 
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Figure 4. TASlc is a master tasiRNA-generating unit regulated by feedback loops. (A) Genetic requirements for the biogenesis of TASl/2-denved 
siRNAs. The panel shows blot hybridization analysis (15% PAGE) of total RNA from CaMV-infected (+) or mock-inoculated (— ) wild-type 
Arabidopsis (Col-0), single (dcl2-5, dcl3-l and dcl4-2) and double (dcl2/3, dcl2/4 and dctf/4) DCL null mutants, AGO mutants (agol-25 and 
ago2-l) and TASlc siRNA-deficient mutant (taslc-1) plants. The cropped scans represent blots successively hybridized with probes specific for 
TASlc-, TASlb-, TASla/b-, TASla/b/c-, TAS2- and TAS3a-derived siRNAs and miR173. EtBr-stained rRNA is shown as a loading control. The 
relative levels of TAS siRNA accumulation in a mutant line versus Col-0 were calculated using the miR173 level as an internal control for each 
sample; for each blot, siRNA levels in Col-0 were set to 1. (B) Blot hybridization analysis of 7MS/c-derived dsRNA-D6 and dsRNA-DlO. The 
cropped scans are from the blot shown in panel (A) (the second lane from left). (C) Blot hybridization analysis of TASlc-, TASlb-, TASla-, TAS2- 
derived dsRNAs in CaMV-infected Col-0 and taslc-1 plants using 5% PAGE. The total RNA samples analyzed are the same as in panel A (the two 
from right). EtBr-stained RNA is shown as a loading control. Positions of dsRNAs derived from TASlc, TASlb. TASla and TAS2 are indicated by 
arrows (same as in Figure IB). TASlc 3'D6(— ) tasiRNA-dependent dsRNAs are highlighted in pink, TASlc 3'D10(— )-dependent dsRNA is in green 
and tasiRNA-independent longest dsRNAs are in blue. (D) Model for TASlc tasiRNA biogenesis and feedback regulation. 



(4)] showed that this 21-nt RNA is one of the most 
abundant siRNAs derived from TASlc. It is shifted 
forward by 1 nt with respect to the 21-nt phase set by 
miR173 cleavage, because a 22-nt siRNA/siRNA* 
duplex is generated at the fourth or the fifth processing 
cycle (Supplementary Figure S3A). The 1-nt phase 
forward drift was also found to occur in TAS3a and 
may result from misprocessing by DCL4 (4). 

Consistent with the deep-sequencing data, our blot hy- 
bridization showed that TASlc 3'D6(— ) appears as one 
21-nt class, while 3'D5(— ) accumulates as minor (21 nt) 
and major (22 nt) classes (Figure 4A). Genetic analysis 
revealed that the biogenesis of both 3'D6(— ) and 
3'D5(— ) requires DCL4, with no substantial contribution 
of DCL2, the dicer normally involved in generating 22-nt 
RNAs (Figure 4A). 



The 3'D6(-) is a 21-nt 5'U-RNA which would predict its 
association with AGOl (15). Indeed, our IP experiment 
using native AGOl -antibodies revealed that AGOl 
protein binds this siRNA as tightly as miR173 (Figure 5A). 
We conclude that both miR173 (18) and TASlc 3'D6(-) 
tasiRNA guide AGOl -mediated cleavage of TAS2 and 
TASlajb pri-RNAs at their respective target sites. 

We also identified additional target sites for TASlc 
3'D6(-) on TAS2 and TASla pri-RNAs, each located 
between the two mapped cleavage sites (Figure 3). 

TASlc 3'D6(— ) directs in trans cleavage of TASlc 
pri-RNA downstream of the origin site, which leads to its 
shorter dsRNA precursor 

TASlc 3'D6(— ) has the potential to cleave its own precur- 
sor (TASlc pri-RNA) in cis, i.e. at the site of its origin 
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Figure 5. Loading of TASlc siRNAs into AGOl and AG02 complexes that mediate in trans and in cis cleavages. (A) The upper panel shows blot 
hybridization analysis of total sRNAs (input) and sRNAs associated with AGOl protein in mock-inoculated (— ) or CaMV-infected (+) plants 
following IP with AGOl -antibodies (IP:AG01). The scans represent the same blot successively hybridized with probes specific for TASlc siRNAs, 
miR173 and, as a negative control, the transposon SIMPLEHAT2-deri\ed 24-nt siRNA. The lower panel shows Western blot analysis of AGOl 
protein accumulation in the input and the IP:AG01 fractions using AGOl-antibody ((a AGOl). Underneath is the same blot stained with antibodies 
specific for CaMV TAV protein (@TAV); (B) The in cis cleavage site for TASlc-deri\ed siRNA 3'D10(— ) at TASlc pri-RNA (the numbering is from 
the 5' cap-site). Thick and thin arrows indicate the major and minor cleavage products, respectively, determined by cRT-PCR; the number of clones 
is given when more than one clone had the same 3'-terminus. The sequence of 3'D10(— ) siRNA guiding the cleavage is shown under the pri-RNA 
sequence. The 5'-termini of predominant antisense RNAs of 1c dsRNA-DIO are shown, with corresponding numbers of specific cRT-PCR products/ 
total products indicated; (C) Blot hybridization analysis of TASlc-derived dsRNAs in CaMV-infected Col-0 and ago2-l plants using 5% PAGE. The 
blot was hybridized with TASlc 3'D6(— ) siRNA-specific probe. EtBr-stained RNA is shown as a loading control. Positions of 7^4S7r-derived 
dsRNAs and dsRNA-DIO are indicated by arrows; (D) sRNA sequencing reads for TASlc-deiived siRNAs and miR173 in sRNA libraries prepared 
from AGO:IP complexes [(15); GEO accession: GSE10036]. 



(Figure 3), which together with miR173-cleavage would 
give rise to a 115-nt fragment. However, we did not 
detect any abundant sense or antisense RNA of this size 
with TASlc probes (Figures 1 and 4B). Instead, we 
identified another complementary site for TASlc 3'D6(— ) 
in a downstream position of TASlc pri-RNA, and con- 
firmed by cRT-PCR that it is cleaved at the expected 
position (Figure 3 and Supplementary Figure S2C). The 
255-nt fragment (U378 to G632) excised from TASlc 
pri-RNA by miR173 and TASlc 3'D6(-) is converted to 



dsRNA with 1 or 2-nt 3'-overhangs on both ends, since the 
termini of the corresponding antisense RNA were mapped 
to positions C631 or U630 (5' end) and A378 (3'-end); two 
non-template nucleotides were added in most cases (Figure 
3 and Supplementary Figure S2C). A sense/antisense RNA 
pair of the corresponding size is the second most abundant 
TASlc pair detected by blot hybridization (Figures 1 and 
4B; Supplementary Figure S2C; designated 1c 
dsRNA-D6). This finding is indicative of a feedback loop 
regulating TASlc 3'D6(— ) biogenesis (Figure 4D). 
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The most abundant precursor of TAS1 3'D6(— ) is 
generated via in cis cleavage 

The most abundant sense/antisense RNA pair derived 
from TASlc is shorter than the earlier described lc 
dsRNA-D6 (Figures 1 and 4B; designated lc 
dsRNA-DIO). cRT-PCR (Supplementary Figure S2C) 
revealed that it is generated by a two-cleavage mechanism 
involving miR173 and TASlc-derived siRNA 3'D10(-) 
which appears to cleave TASlc pri-RNA in cis at the 
sites of its origin (Figure 5B). The sense strand of the 
resulting dsRNA starts with U379 and terminates pre- 
dominantly with U577. The antisense strand 5'-ends map 
predominantly to U576, U575 and A574, thus resulting 
in 1- to 3-nt 3'-overhang (Figure 5B). The corresponding 
3'-ends map to A379 followed with a few non-template 
nucleotides (Supplementary Figure S2C). 

The 3'D10(— ) itself can be generated from both TASlc 
longest dsRNA and lc dsRNA-D6 (Figure 4C) by DCL4 
processing at the tenth cycle from the miR173 end. On 
blots, the 3'D10(— ) appears as 21-nt and, more 
abundant, 22-nt classes (Figure 4A). The mapped 
cleavage site on TASlc pri-RNA would predict that 
both size-classes should have 5'-terminal A587 to direct 
an AGO-mediated cleavage at the mapped position 
(Figure 5B). 

Genetic analysis showed that both 21-nt and 22-nt 
classes of TASlc 3'D10(— ) require DCL4 for their biogen- 
esis. The 21-nt class is selectively reduced in the null 
mutant ago2-l (15), whereas the 22-nt class is selectively 
reduced in the hypomorphic mutant agol-25 (26) 
(Figure 4A). This suggests that TASlc 3'D10(-) exists 
as 21-nt 5'A587-RNA bound to AG02 and 22-nt 
5'U588-RNA bound to AGOl. Our AGOl-IP and 
other AGO-IP studies support this conclusion 
(Figure 5 A and D). 

Taken together, our findings suggest that the 21-nt 5'A- 
form of TASlc 3'D10(-) may guide AG02-mediated 
cleavage of TASlc pri-RNA in cis. AG02 is known to 
be associated with 21-nt 5'A-RNAs (15) but a slicer 
activity of AG02 was not demonstrated so far. 

In ago2-l mutant, apart from the reduced accumulation 
of 21-nt 3'D10(— ), accumulation of other 7^,S7c-derived 
tasiRNAs including 3'D6(— ) is not drastically affected 
(Figure 4A). This can be explained by increased (compen- 
satory) production of tasiRNAs from lc dsRNA-D6 and 
also from TASlc longest dsRNA. Indeed, accumulation 
of lc dsRNA-D6 is increased ~2- to 3-fold in ago2-l 
compared to wild-type plants, concomitant with ~2- to 
3-fold decrease in accumulation of 7c-dsRNA-D10; the 
longest dsRNA remains at a similar low level (Figure 5C 
and Supplementary Figure S4). Thus, in the absence of 
AG02, lc dsRNA-D6 becomes the most abundant pre- 
cursor of 7L4,S7c-derived tasiRNAs. 

Surprisingly, in dc/2/3/4 plants that do not produce any 
tasiRNA, in addition to the longest TASlc dsRNA, one 
shorter dsRNA accumulates to high levels (Figure 1). This 
dsRNA has 3'-overhangs of one to several nucleotides and 
is generated by a mechanism involving the miR173 
cleavage and a precise cleavage that occurs predom- 
inantly 7-nt upstream of the 3'D6(— ) cleavage site 



(Supplementary Figure S2C). It remains to be investigated 
which DCL1 -dependent sRNA directs the second cleavage 
in this case. 

The majority of TAS1/TAS2 siRNAs are produced from 
the region between the cleavage sites 

Our inspection of the deep-sequencing data (4) revealed 
that for each miR173-dependent TAS gene, most of 
unique sRNAs and the overwhelming majority of sRNA 
reads originate from the region between miR173 and 
TASlc 3'D6(— ) cleavage sites (Supplementary Figure 
S3B). Thus, the short dsRNAs (dsRNA-D6) produced 
by the two-cleavage mechanism are likely the major pre- 
cursors of TAS1/TAS2 tasiRNAs. The region of TASla/ 
b/c and TAS2 genes downstream of the TASlc 3'D6(— ) 
cleavage site also generates sRNAs, albeit at much lower 
levels (Supplementary Figure S3B). These sRNAs are 
likely processed from the longest, poly(A)/poly(U) 
duplex-containing dsRNAs. 

TAS2- and TASlalb/c-demed short dsRNAs are 
processed from both ends 

The dsRNAs derived from TAS2, TASla, TASlb and 
TASlc can potentially be processed from both ends. If 
this is true, then in addition to miR173-phased 
tasiRNAs, siRNAs in phase with the TASlc 3'D6(-) 
cleavage site should also accumulate. Inspection of the 
deep-sequencing data (4) revealed that each of the four 
miR173-dependent TAS loci spawns 21-nt RNAs 
upstream of and in phase with the 3'D6(— ) cleavage site 
(Supplementary Figure S3A; highlighted in pink). Also, 
more abundant 21-nt RNAs are produced in a 1-nt- 
forward-drifted phase caused by generation of 22-nt 
RNAs at the first {TASla and TASlb) or the seventh pro- 
cessing cycle (TAS2) (Supplementary Figure S3A; purple). 
Thus, similar to miR390-directed cleavage for TAS3 
family pri-RNAs (12), the 3'D6(-) cleavage of TASla/b/ 
c and TAS2 pri-RNAs sets a register for production of 
phased siRNAs upstream of the cleavage site. This 
second register is shifted with respect of the miR173 
register either forward (by 2nt in TAS2 and TASlc and 
3nt in TASla) or backward (by 3 nt in TASlb). 
Interestingly, an expected cleavage of TASla and TAS2 
pri-RNAs at the additional target sites for TASlc 3'D6(— ) 
(Figure 3) also appears to initiate processing of phased 
siRNAs from the upstream region (Supplementary 
Figure S3A; cyan and blue). 

Likewise, TASlc 3'D10(-) cleavage of TASlc pri-RNA 
sets a register for 21-nt RNAs from the upstream region 
(Supplementary Figure S3A; green). Notably, the 22-nt 
TASlc D10(— ) species itself can be produced by 
in-phase processing of lc dsRNA-D6 from either end. 

Using a miRNA target prediction software [http:// 
bioinfo3.noble.org/psRNATarget/; (27)] we identified po- 
tential target genes for several siRNAs encoded in the regis- 
ters set by the 3'D6(-) on TASla, TASlb and TAS2 and 
the 3'D10(-) on TASlc (Supplementary Table SI). 

Taken together, we conclude that the short dsRNAs 
derived from TASla, TASlb, TASlc and TAS2 are pro- 
cessed by DCL4 not only from the miR173 end but also 



6250 Nucleic Acids Research, 2012, Vol. 40, No. 13 



from the opposite end defined by TASlc 3'D6(— ) cleavage 
and, in the case of TASlc, also by TASlc 3'D10(-) 
cleavage. The processing from both ends results in produc- 
tion of siRNAs in two different registers, thus increasing 
the coding capacity of these tasiRNA-generating loci. 

It should be noted that the deep-sequencing data used 
for the above analyses of tasiRNAs were obtained using 
non-infected A. thaliana plants (4). Thus, the short 
dsRNA precursors of tasiRNAs identified here in 
CaMV-infected plants are likely present and processed 
from both ends in non-infected plants. 

miR173-dependent tasiRNA production from TASla/b 
and TAS2 genes is not drastically affected by TASlc 
gene knockdown 

Our results so far demonstrated that 7L4S7c-derived 
tasiRNA 3'D6(-) can target in trans TASla/b and TAS2 
pri-RNAs to initiate the biogenesis of abundant dsRNAs 
that can be processed by DCL4 from both ends. To address 
the question whether this second cleavage is required for 
efficient production of TASla/b and TAS2 tasiRNAs in the 
miR173 register, we searched A. thaliana mutant collec- 
tions for TASlc gene mutants. By genotyping combined 
with RNA blot hybridization, we identified a homozygous 
mutant line (designated taslc-1) with T-DNA insertion 
between the two major poly(A) sites of TASlc 
(Supplementary Figure S2C). In taslc-1, accumulation of 
TASlc siRNAs including 3'D6(-) and 3'D10(-) was 
reduced to ~40% of the wild-type level in rosette leaves 
and to ~10% in inflorescence tissues (Figure 4A). This 
knockdown of TASlc gene expression strongly reduced ac- 
cumulation of the 3'D6(— )-dependent short dsRNAs 
derived from TASla, TASlh and TAS2 (Figure 4C; see 
'dsRNA-06 1 species highlighted in pink), but did not 
affect the accumulation the 3'D6(— )-independent, longest 
dsRNAs derived from TASla and TASlh (Figure 4C; high- 
lighted in blue) as well as less abundant, shorter RNAs 
derived from TAS2 (Figure 4C; indicated by four con- 
nected arrows). The knockdown had only a slight impact 
on accumulation of TASla/b and TAS2 tasiRNAs encoded 
in the miR173 register, which was reduced only to 70-80% 
of the wild-type levels in both tissues (Figure 4A). We 
conclude that TASlc gene expression does not substan- 
tially affect the biogenesis of TASla, TASlb and TAS2- 
derived siRNAs encoded in the miR173 register. Our 
findings are consistent with the notion that tasiRNA pro- 
duction from TASla/b and TAS2 loci can proceed through 
a single miR173-directed cleavage event (18,19), which 
leads to the production of the longest, poly(A)/poly(U) 
duplex-containing dsRNAs. 

TASlc 3'D6(— ) tasiRNA-directed cleavage regulates 
TASlc 3D10(-) siRNA biogenesis and precludes TAS2 
siRNA production downstream of the target site 

To further investigate a role of the second cleavage of TAS 
pri-RNAs, we used a previously established transient ex- 
pression assay using N. benthamiana leaves which recap- 
itulates miR173-dependent tasiRNA biogenesis from the 
TAS1 and TAS2 transcripts driven by the CaMV 35S 
promoter (10). Confirming the previous results, robust 



production of siR255 from 35S:TASlc and siR1940 
from 35S:TAS2 could be achieved only by co-expression 
of miR173 from a 35S:miR173 construct. We then 
co-expressed the TAS2 construct with 35S:TASlc or its 
derivative in which the seed sequence of the 3'D6(— ) 
tasiRNA had been mutated to prevent the second 
cleavage of TAS2 pri-RNA. Co-expression of 3'D6(— ) 
species or its seed-mutant form only slightly affected 
miR173-dependent production of siR1940 from 7^5*2 
(Figure 6), supporting the above conclusions. However, 
siRNA production from the region downstream of the 
3'D6(— ) cleavage site was nearly abolished in the 
presence of 3'D6(— ) species but not its mutant form (see 
'TAS2 D6 down 1 in Figure 6). Thus, the second cleavage 
of TAS2 pri-RNA restricts siRNA production to the 
region between the miR173 and 3'D6(— ) cleavage sites. 
This is consistent with our above analysis of the sRNA 
deep-sequencing data in A. thaliana. Interestingly, the 
region downstream of the 3'D6(— ) cleavage site generated 
siRNAs of three size-classes and their production was 
detected both in the presence and absence of miR173. 

For TASlc itself, the mutation of the 3'D6(— ) seed 
sequence did not alter miR173-dependent accumulation 
of siR255. However, it reduced accumulation of 22-nt 
3'D10(— ) and increased accumulation of 21-nt 3'D10(— ) 
species, thus resembling the effect observed in Arabidopsis 
agol-25 plants. Mutation of the 3'D6(— ) target site (TS) in 
35S:TASlc, in such a way that the mutant TS sequence 
could be cleaved only by the D6 seed-mutant siRNA but 
not the wild-type siRNA, led to a similar effect on accu- 
mulation of the 3'D10(— ) size-classes. A combination of 
the two mutations restored the wild-type ratio of 3'D10(— ) 
size-classes. Furthermore, an A-to-G mutation of the 5'- 
terminal nucleotide of the 21-nt 3'D10(-) in 35S:TASlc 
increased the ratio of the 22-nt class to the 21-nt class, 
which is reminiscent of the effect observed in 
Arabidopsis ago2-l plants. Taken together, these results 
confirm feed-back regulation of TASlc siRNA biogenesis 
by the 3'D6(— )-directed in trans cleavage and the 
3'D10(— )-directed in cis cleavage of TASlc pri-RNA. 

DISCUSSION 

Structures of RDR6-dependent dsRNA precursors of 
tasiRNAs 

In this work we demonstrate that RDR6 converts the 
entire, polyadenylated, 3'-products of miR173-cleaved 
TAS1/TAS2 pri-RNAs to dsRNA. However, the resulting 
longest dsRNAs accumulate at much lower levels than the 
shorter dsRNAs produced from all the four loci by a 
two-cleavage mechanism. In fact, 7^4S2-derived longest, 
poly(A)/poly(U) duplex-containing dsRNA could be 
detected only in dcl2/3/4 plants lacking DCL2, DCL3 
and DCL4 activities, while TA SI /oc/-derived longest 
dsRNAs accumulated at much higher levels in dcl2/3/4 
than in wild-type plants. Thus, the polyadenylated 
3'-fragment of miR173-cleaved TAS pri-RNAs is rapidly 
truncated by the second cleavage directed by TASlc- 
derived tasiRNA 3'D6(— ) which generates a shorter, 
poly(A)-less template for RDR6. We cannot exclude, 
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Figure 6. TASlc 3'D6(— ) tasiRNA-directed cleavage regulates TASlc 3'D10(— ) siRNA biogenesis and precludes TAS2 siRNA production down- 
stream of the target site. (A) The wild-type and mutant targets sites (TS) of the wild-type TASlc pri-RNA and its TS mutant version, respectively. 
Mutagenized positions are in lower case and highlighted in cyan. Position of the cleavage directed by the wild-type TASlc 3'D6(— ) tasiRNA or its 
seed-mutant version is indicated by red arrow. (B) Reconstruction of TASlc and TAS2 tasiRNA biogenesis in a transient expression assay using 
N. benthamiana. 35S:TASlc constructs with mutagenized 3'D6(— ), 3'D10(— ) siRNAs, the 3'D6(— ) target site sequences, or both 3'D6(— ) and its 
target site sequences were expressed or coexpressed with 35S:miR173 and also with 35S:TAS2 constructs as indicated above the blot panels. 
The siRNAs detected by successive hybridization of the blots are shown to the right of each panel. 



however, that the longest dsRNA is more efficiently pro- 
cessed by DCL4 than the shorter dsRNA. 

Our cRT-PCR analysis also suggests that the poly(A) 
tail is shortened and, in some cases, fully removed prior to 
RDR6 binding to the long template. The shortened 
poly(A) tails are often tagged by short UA-rich sequences. 
The RNA templates with both tagged and non-tagged 
poly(A) tails of up to 37 nt are copied by RDR6 to 
dsRNAs. This finding is supported by the in vitro studies 



of an affinity-purified A. thaliana RDR6 protein showing 
that the RDR6 activity is not inhibited by addition of a 
40-nt poly(A) tail to the template (28). 

The structures of more abundant, short TAS1/TAS2 
dsRNAs identified here reveal that RDR6 initiates com- 
plementary RNA synthesis preferentially from the third 
nucleotide of the template's 3'-end, thus creating a 2-nt 
3'-overhang on the resulting dsRNA. Less frequently, 
the synthesis is initiated from the second or the fourth 
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nucleotide, resulting in 1-nt or 3-nt 3'-overhangs, 
respectively. In our parallel work, precise mapping of 
the TAS3a-devived dsRNA revealed that RDR6 initiates 
complementary RNA synthesis predominantly at the third 
nucleotide of the miR390-cleaved TAS3a pri-RNA (17). 
This confirms our findings for miR173-dependent TASlj 
TAS2 dsRNAs. Furthermore, we demonstrate here that 
RDR6-mediated synthesis terminates at the 5'-terminal 
nucleotide of the template and frequently one or two 
non-template nucleotides are added to the 3'-end of the 
complementary RNA. Similar biochemical properties 
have been reported for an RDR activity isolated from 
tomato: this RDR copied an artificial single-stranded 
RNA template from the second or the third nucleotide 
of the template's 3'-terminus carrying the hydroxyl 
group and added one non-template nucleotide to 50% of 
the complementary run-off transcripts (25). Thus, produc- 
tion of dsRNAs with 1-nt or 2-nt 3' overhangs at both 
termini appears to be a common property of plant 
RDRs. In vitro studies have shown that the dsRNA with 
a 2-nt 3'-overhang is an ideal substrate of human dicer 
(29). It remains to be investigated whether or not a short 
3' overhang is absolutely required for DCL4. Notably, a 
substantial fraction of miR173-dependent dsRNAs 
derived from TASla, TASlb, TASlc or TAS2 does not 
possess any overhang at the miR173 cleavage end. 

Our data also support the in vitro results showing that 
Arabidopsis RDR6 does not require a primer (28). 

Role of the second cleavage of pri-RNA in tasiRNA 
production from TAS1/TAS2 genes 

We demonstrate here that miR173-initiated biogenesis of 
tasiRNAs from Arabidopsis TAS2 and TAS1 genes 
involves a second cleavage of TAS pri-RNA downstream 
of the miR173 site. In all cases, the short segment of 
pri-RNA between the two-cleavage sites is converted by 
RDR6 to dsRNA. The resulting dsRNA spans the region 
giving rise to the majority of siRNAs derived from each of 
the four loci. Hence, the two-cleavage mechanism gener- 
ates the major dsRNA precursors of tasiRNAs. However, 
our results also demonstrate that the second cleavage is 
not required for production of tasiRNAs in the miR173 
register. This is consistent with the previous studies 
showing that a single miR173-dependent cleavage is suffi- 
cient to initiate tasiRNA biogenesis (18,19). What could 
be the role of the second cleavage? Our results in 
N. benthamiana indicate that the second cleavage 
directed by TASlc 3'D6(-) tasiRNA strongly reduces 
TAS2 siRNA production from the region downstream 
of the cleavage site. Thus, the two-cleavage mechanism 
appears to restrict siRNA production to the region 
between the miR173 and 3'D6(— ) cleavage sites. 
Additionally, identification of 2 1-nt phased siRNAs in 
the register set by TASlc 3'D6(— ) cleavage in the sRNA 
deep-sequencing database shows that the second cleavage 
could also be to increase the tasiRNA coding capacity of 
the TAS1/TAS2 genes. 

Our findings extend the two-hit trigger model originally 
proposed for miR390-dependent TAS3 genes (12) to 
miR173-dependent TAS1/TAS2 genes. But unlike 



TAS3a/b/c pri-RNAs targeted at two sites by a single 
miRNA (miR390), TASla/b/c and TAS2 pri-RNAs are 
targeted by miR173 only at one site, while the second 
hit is guided by 7MS7c-derived tasiRNA that cleaves not 
only TASla, TASlb and TAS2 pri-RNAs but also its own 
precursor (TASlc pri-RNA). Previous studies have 
identified several Arabidopsis genes that have target sites 
for two or more different miRNAs and/or tasiRNAs; 
most of these genes generate secondary siRNAs (4,12). 
We demonstrate here that the miRNA- and the 
tasiRNA-directed cleavages at one target RNA can 
recruit RDR6 activity that converts the RNA segment 
between the two-cleavage sites to a dsRNA substrate for 
DCL4. 

Feedback regulation of TASlc gene 

We propose a feedback regulation model in which TASlc 
gene codes for two siRNAs, 3'D6(— ) and 3'D10(— ), which 
control tasiRNA production from this gene (Figure 4D). 
These siRNAs are encoded in antisense orientation at the 
processing cycles 6 and 10, in phase with the miR173 
register, and can initially be generated from the poly(A)/ 
poly(U) duplex-containing longest dsRNA, produced by 
RDR6 following a single miR173-cleavage. Once 
generated, each of these siRNAs together with miR173 
directs cleavage of TASlc pri-RNA at the respective 
target site to initiate RDR6-catalyzed synthesis of two 
distinct dsRNAs. The most abundant of the two (7c 
dsRNA-DIO) is generated by the 3'D10(-) cleavage in 
cis and can potentially serve as a precursor of 3'D6(— ) 
siRNA but not of 3'D10(-) itself. Less-abundant 
dsRNA (lc dsRNA-D6) generated by the 3'D6(-) 
cleavage in trans can potentially spawn both species 
(Figure 4D). 

Several lines of evidence support our feedback model 
and allow speculation on potentially positive or negative 
functions of the two feedback loops. First, we have 
detected by blot hybridization and sequenced the three 
7/lS7c-derived dsRNAs, namely, the single 
miR173-dependent longest dsRNA (7c dsRNA-long; 
Figures IB and 2), the miR173- and TASlc 3'D6(-) 
tasiRNA-dependent dsRNA (7c dsRNA-D6; Figures IB 
and 3) and the miR173- and TASlc 3'D10(-) 
siRNA-dependent dsRNA (7c dsRNA-DIO; Figures IB, 
3 and 5). The terminal sequences of the dsRNA template 
strands indicate the conserved, AGO-mediated cleavages 
of TASlc pri-RNA directed by miR173 on the upstream 
site and by the 3'D6(-) or 3'D10(-) siRNA at the down- 
stream site. The relative accumulation levels of the three 
dsRNAs would suggest that DCL4 processing of the two 
shorter dsRNAs might generate the majority of TASlc- 
derived siRNAs. However, we cannot exclude that the 
longest dsRNA is processed more efficiently than the 
shorter ones. Second, we provide genetic evidence that 
reduced accumulation of 21-nt 3'D10(— ) siRNA in 
ago2-l mutant plants correlates with reduced accumula- 
tion of 7c dsRNA-DIO and increased accumulation of 7c 
dsRNA-D6 (Figure 5). Thus, the 3'D10(-)-directed 
cleavage required for 7c dsRNA-DIO biogenesis appears 
to interfere with 7 c dsRNA-D6 biogenesis, possibly owing 
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to competition between 3'D10(-) and 3'D6(-) siRNAs 
for the same target TASlc pri-RNA. This suggests that 
the in cis cleavage directed by 3'D10(— ) may be involved 
in negative feedback regulation of TASlc siRNA biogen- 
esis. Third, our reconstitution of TASlc siRNA produc- 
tion in the N. benthamiana transient assay shows that the 
3'D6(-) tasiRNA-directed cleavage of TASlc pri-RNA 
regulates the biogenesis of 3'D10(— ) siRNA size-classes 
by decreasing relative accumulation of the 
cleavage-competent, 21-nt species (Figure 6), thus suggest- 
ing a positive function of this cleavage event. Further 
experiments are needed to clarify the roles of the two 
feedback loops in regulation of TASlc gene expression. 

We established that, similar to miR173, TASlc 3'D6(-) 
tasiRNA is tightly associated with AGOl protein to 
exert its function in TAS1/TAS2 siRNA biogenesis. 
Interestingly, this 21-nt 5'U-RNA was also 
co-immunoprecipitated with AG02 which preferentially 
binds 21-nt 5'A-RNAs and with AG04 which preferen- 
tially binds 24-nt 5'A-RNAs [(15); Figure 5D]. We specu- 
late that this apparent competition for the 3'D6(— ) species 
may play a role in regulation of TASlc expression by 
reducing the levels of functional AG01::3'D6(— ) 
complexes. Additionally, AG04::3'D6(— ) complex may 
function in the nucleus to downregulate transcription of 
TASlc pri-RNA. 

AG02 was implicated in TAS3 siRNA biogenesis and 
in antiviral defense since it was found to be associated 
with miR390 (15) and virus-derived siRNAs (16), respect- 
ively. Recently, the antiviral role of AG02 was confirmed 
for some RNA viruses that encode suppressors specifically 
targeting AGOl (30). However, so far there was no 
evidence for AG02 catalytic activity. Our results suggest 
that this protein binds 7^457c-derived 21-nt siRNA 
3'D10(— ) to mediate in cis cleavage of TASlc pri-RNA. 
Further experiments are needed to confirm that AG02 is a 
functional slicer and that its sheer activity is required for 
AG02 action in TASlc tasiRNA biogenesis and antiviral 
defense. 

TAS1/TAS2 gene regulation network 

A sRNA interaction network connecting miR173- 
dependent TAS2, TASla and TASlc genes was recently 
predicted by a bioinformatic analysis of sRNA popula- 
tions sequenced from A. thaliana (6). Our data confirms 
the existence of this network as well as refines and extends 
the prediction by demonstrating that 7^4,S7c-derived 
tasiRNA can target all the four miR173-dependent TAS 
loci. Few unmapped RNAs that we detected by blot hy- 
bridization (Figure 1) may arise from additional targeting 
events predicted by MacLean et al. (6). 

The existence of this complex network with a central 
role for TASlc is not surprising, given the fact that the 
four 7^45 loci are paralogous: TAS2 and TASlc arose by 
a direct duplication, while TASla and TASlb arose from 
TASlc after it had diverged from TAS2 (11). Interestingly, 
TAS2 and TAS1 tasiRNAs including TASlc 3'D6(-) and 
3'D10(— ) described here target several pentatricopeptide 
repeat family genes giving rise to secondary siRNAs with 
the potential to target in trans other genes within and 



outside of this multigene family (4). This creates another 
complex interaction network subordinate to the miR173 
network. Our findings for TASlc would also suggest that 
feedback loops in the hubs of sRNA networks can 
fine-tune the levels of a master siRNA, such as TASlc 
3'D6(— ), required for in trans regulation of other genes. 
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